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An explicit expression is derived for the scattering function of a self-avoiding polymer chain in a 
d-dimensional space. The effect of strength of segment interactions on the shape of the scattering 
function and the radius of gyration of the chain is studied numerically. Good agreement is demon- 
strated between experimental data on dilute solutions of several polymers and results of numerical 
simulation. 
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I. INTRODUCTION 

This paper is concerned with the scattering function 
(form factor) of a self-repellent polymer chain. The ef- 
fect of excluded-volume interactions between segments 
on the characteristic size of a "real" chain has attracted 
substantial attention from the beginning of the 1950s 
0,1^1^. Although a noticeable progress has been reached 
in this area (the results of investigation are summa- 
rized in monographs 0, 0, IE S Hi) i the effect of self- 
avoiding interactions on the statistics of macromolecules 
has remained a subject of debate in the past decade, see 
H ES in 111 III Hi El El IS, to mention a few. 
This fact may be ascribed to the importance of excluded- 
volume interactions for the description of such phenom- 
ena as (i) protein folding and denaturation J^, JSj , (ii) 
unzipping of DNA |0| , (iii) mechanical response of indi- 
vidual macromolecules , and (iv) transport of chains 
through narrow pores jig . A correct account for self- 
repellent interactions (i) provides a consistent descrip- 
tion of the effect of molecular w eigh t on intrinsic viscos- 
ity of dilute polymer solutions E3|, (ii) leads to an ad- 
equate characterization of time-dependent shear stresses 
im ll^ and transport coefficients in polymer fluids |0| , 
(iii) predicts reinforcement of interfaces between incom- 
patible solvents by random copolymers |^, etc. 

The influence of self-avoiding interactions on the scat- 
tering intensity of a chain in light, X-ray, and neu- 
tron scattering experiments has been investigated in 
10., J3, J4, 25, Jfi], where phenomenological 

equations have been developed for the structure func- 
tion based on simplified theoretical models and results 
of Monte Carlo simulation. Although these relations 
improve our knowledge of the statistics of self-avoiding 
chains, their applicability remains questionable for the 
lack of direct comparison of their predictions with those 
grounded on an exact expression for the form factor. The 
objective of this study is to derive analytical formulas 
for the scattering function of a self-avoiding chain in a 
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d-dimensional space and its radius of gyration and to 
analyze the influence of strength of excluded-volume in- 
teractions on these quantities. 

Two approaches are conventionally employed to model 
a self-repellent flexible chain. According to the first, a 
chain is treated as a set of freely-jointed rigid segments 
linked in sequel, whereas self-avoiding interactions do not 
ermit each pair of segments to occupy the same place 
l|. This model implies that the distribution function of 
end-to-end vectors for the chain coincides with that for a 
self-avoiding random walk on a lattice |2l 13 ■ Although 
the latter function is unknown, it is believed that the 
leading term in the expression for the distribution func- 
tion of a long chain is described by a stretched exponen- 
tial function multiplied by a polynomial of the end-to-end 
distance pllsot. 

The other approach (the Edwards model lU) treats 
an arbitrary configuration of a chain as a curve with a 
Hamiltonian ascribed to it. The Hamiltonian equals the 
sum of the conventional Hamiltonian for a Gaussian chain 
and the excluded-volume functional, where the energy 
of segment interactions is approximated by the delta- 
function. This method does not exclude intersections 
between segments, but assumes their number to be re- 
duced (compared to that for a Gaussian chain) in some 
proportion to the strength of self-repellent interactions. 
The distribution function of end-to-end vectors for a self- 
avoiding chain is determined by the conventional path 
integral (Green's function), and the scattering function 
(form factor) is given by the classical formula |33| , where 
averaging is performed with the help of the Green func- 
tion. 

Despite an apparent simplicity of the Edwards model, 
no exact formulas have yet been derived within this 
approach for the radial distribution function and the 
form factor. This may be explained by the fact that 
the Hamiltonian containing the excluded-volume func- 
tional becomes non-quadratic, and the path integral can- 
not be computed by using conventional approaches |33j| . 
Among non-rigorous techniques proposed for the treat- 
ment of t his p roblem, it is worth mentionin g pe rturba- 
tive HIM Hi Hi and renormafization-group |37L Isl Isgf 
methods. 
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The purpose of the present paper is to derive an ex- 
act formula for the form factor of a self-avoiding chain 
in a d-dimensional space based of the Edwards model 
with a slight modification of the excluded-volume func- 
tional. This modification is caused by the observation 
that the average (over configurations) value of the Ed- 
wards Hamiltonian, as well as all moments of the number 
of segment intersections become infinite at d > 1. The 
divergence reflects an increase in the number of intersec- 
tions between different segments of the chain when the 
distance between them (measured along the curve that 
describes an arbitrary configuration) approaches zero. 
Based on a similarity between configurations of a Gaus- 
sian chain and trajectories of a standard Brownian mo- 
tion, the latter may be explained by the well-known prop- 
erties of the local time of intersection for a d-dimensional 
Brownian path ■ 

We concentrate on the derivation of an explicit ex- 
pression for the scattering function. The interest to this 
quantity is explained by the fact that a number of experi- 
mental data on small-angle X-ray and neutron scattering 
by macromolecules is available, which allow the strength 
of excluded-volume interactions to be found by fitting 
observations. An analytical formula for the distribution 
function of end-to-end vectors will be developed in a sub- 
sequent paper. The exposition is organized as follows. 
The average number of self-intersections for a flexible 
chain is calculated in Section 2. The Edwards model with 
a regularized Hamiltonian is introduced in Section 3. A 
series expansion for the scattering function is developed 
in Section 4. Borel's summation of the series is performed 
in Section 5, where a closed- form expression is derived for 
the Laplace transform of this function. Explicit expres- 
sions for the scattering function of a self-repellent chain 
are obtained in Section 6 for the one-dimensional case 
and in Section 7 for the three-dimensional case. Two 
asymptotic approaches are suggested for the simpliflca- 
tion of the form factor for a self-avoiding chain. The 
approximate expressions are tested in Section 8, where 
results of numerical simulation are compared with exper- 
imental data. Some concluding remarks arc formulated 
in Section 9. 



II. THE LOCAL TIME OF 
SELF-INTERSECTIONS 

We consider a standard Gaussian chain in a d- 
dimensional space. One end of the chain is flxed at the 
origin, and the other end is located at a point with ra- 
dius vector Q. An arbitrary configuration of the chain is 
described by a curve r(r) with t € [0, L], where L stands 
for the contour length. The Hamiltonian of the chain 
reads 



Ho{r) = 




where bo denotes the length of a statistical segment, fee 
is Boltzmann's constant, and T is the absolute tempera- 
ture. 

Fixing a segment of the chain labeled hy t £ (0, L], 
we consider intersections of this segment with other seg- 
ments labeled by t £ (si,S2), where < si < S2 < t. 
The number of intersections (the local time of intersec- 
tion) is defined as 

v{s,,S2;t)^b',-' r <5(r(t) - r(T))dr, (2) 

Jsi 

where 

is the Dirac function, the dot stands for inner product, 
and integration is performed over all vectors k. The 
delta-function in Eq. |5J| counts the "number" of inter- 
sections between the segments labeled by t and r, and 
the integral describes the sum these "numbers" over var- 
ious segments t belonging to the interval (si,S2). The 
pre-factor bo is introduced to make the right-hand side 
of Eq. ^ dimensionless. The dependence of on the 
end-to-end vector Q is suppressed. 

Although Eq. (O is widely used, it should be treated 
with caution for this formula has an unambiguous mean- 
ing when the average (over configurations) number of 
self-intersections (z^(si, S2; i))Q is finite for any point t 
and an arbitrary interval (si,S2). This is not the case, 
however, because all moments of the function ^{si, S2] t) 
diverge as S2 t. This fact appears to be natural, for all 
configurations are taken into account in the calculation 
of the average value of f, including those curves r(T) that 
"oscillate" (and intersect themselves) with high frequen- 
cies. 

Several approaches were developed to avoid the diver- 
gence of the integral in Eq. The simplest method 
|4r| consists in (i) replacement of the delta-function in 
Eq. 1I2J1 with another delta-function that "measures" the 
number of intersections "shifted" to some vector a, (ii) 
determination of the average number of intersections, and 
(iii) calculation of its limit (if it exists) as a ^ 0. This 
technique allows the first moment of the number of in- 
tersections to be found, but does not permit the entire 
set of moments to be determined analytically. 

The renormalization group approach was applied to 
the analysis of intersections of Brownian motion in 
[4^ E3i EJI • An appropriate renormalization (both mul- 
tiplicative and additive) allows (i) explicit expressions for 
the first two moments of i> to be derived, and (ii) conver- 
gence of the "renormalized" process i>{si,S2,t) [treated 
as a function of t for a given interval (si, S2)] to a Brown- 
ian motion with known characteristics to be proved. Al- 
though these methods seem rather promising, their in- 
terpretation in terms of physical parameters of a chain 
remains unclear. 

To avoid complications of the renormalization-group 
approach, we regularize the number of intersections by 
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ascribing to each intersection between segments labeled 
by T and t some weight p that (i) approaches zero when 
the difference t — t vanishes, and (ii) reaches its ultimate 
value of unity when the "distance" between these points 
along the curve r(T) becomes relatively large. This means 
that we replace definition Q with 



^(si,S2;i) 



,d-i 



p{t - T)d{r{t) - r(r))dT. (4) 



As the regularizing function p{t) may be chosen from a 
set of functions satisfying the conditions p(0) = and 
limt^ioo = 1, Eq- l@J provides an additional de- 
gree of freedom in the design of models for self-avoiding 
chains. In what follows, wc set 



p{t) = l-exp(|<|/Lo), 



(5) 



where Lq is the characteristic length of internal inho- 
mogeneity. Roughly speaking, Eqs. Q and Q treat 
intersections between segments located [along the curve 
r(r)] at distances exceeding Lq with the same (unit) 
weight and disregard intersections between segments lo- 
cated [along this curve] closer than the characteristic 
length Lq. 

To find the average number of self-intersections, we 
replace the delta-function in Eq. with its Fourier 
transform (O , and average the obtained relation over all 
configuration of a chain with an end-to-end vector Q, 



{v{si,S2;t))ci 



,d-i 



p{t ~ T)dr 



cxp 



-ik- 



r(i)-r(r)) 



"o 



(27r)" 



Pit 



2d 



A{t,T) -^k• Q 



dk 

Q 

T)dr 



dk. 



where the function A(t, r) describes correlations between 
segments. Calculation of the Gaussian integral results in 



{v{si,S2;t))Q 



1 /dbp 
ba V 27r 



X exp 



pit 



AHt,T) 



dQ^it 



2boL^A{t,T) 



dr. (6) 



Equation ^ describes the average number of self- 
intersections for a chain with an arbitrary correlation 
function A(t, r). For a standard Gaussian chain with 



A(t,r) = (t-r)(l-^), 



we arrive at the formula 

1 rdbp 
bo \ 27r 



(i/(si,S2;i))Q 



X 1 



exp 



Pit 
dQ^it^ 



-r){t- 
r) 



2boL^il~it~T)/L) 



dr. (7) 



The average (over chains with various end-to-end vectors 
Q) number of intersections is given by 



where 



(i^(si,S2;i)) - / (i/(si,S2;i)>QPo(Q)dQ, (8) 



.2TrboLJ ""^V 2boLJ' ""^^ 

is the distribution function of end-to-end vectors for a 
Gaussian chain in a c?-dimensional space. Insertion of 
expressions JTJ and lO into Eq. ^ and calculation of 
the Gaussian integral yield 

bo V 2tt 



(j/(si,s2;0> 



-dr. 



(10) 



(t-r)\ 

Equation (|10|l shows that the average number of intersec- 
tions between a segment labeled by t and segments be- 
longing to the interval (si,S2) approaches infinity when 
■32 t and d > 1, if a regularization function p{t) is not 
used. 

Setting si = and S2 = t in Eq. wc find the 

regularized number of intersections between a segment 
labeled by t and all segments labeled hy t < t. Summa- 
tion of these numbers over t e (0, L] provides us (after 
division by L) with the number of self-intersections per 
unit length 



1 

L Jo 



i^iO,t;t)dt. 



(11) 



Substitution of expression (fTU)! into Eq. ((TT|l results in 
the average number of self-intersections per unit length 
of a standard Gaussian chain in a d-dimensional space 

^ ' boL\27rJ Jo Jo {t~T)i 

where we set s = t — r in the internal integral, changed 
the order of integration, computed the integral over t, 
and set t — s/L in the result. For a Gaussian chain in a 
three-dimensional space, Eqs. lO and (|12() yield 



36n 



2nL 



;^((« + lh(i,»)-7(|«))-2; 



where 



^0 



f z''-iexp(- 

"'0 



z)dz. 



(13) 



When K — > (X), which means that the contour length 
of a chain substantially exceeds the characteristic length 
of regularization, the incomplete gamma-function 7(0, k) 
approaches the gamma- function r(a), and the leading 
term in the expression for the average number of self- 
intersections per unit length reads 



TT V 2Lo 



(14) 
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III. THE EDWARDS MODEL 

According to the Edwards model [sj , the Hamiltonian 
H oi a, self-avoiding chain equals the sum of the Hamil- 
tonian of a standard Gaussian chain Hq and a quantity 
proportional to the number of self-intersections per unit 
length ly*. Equations Q), Q and |[TT|l imply that 



Hir) = ksT 



vU 



2L 



<* L p L 

X I dt p{t - T)S{r{t) - r{T))dT 
h Jo 



, (15) 



where v stands for the strength of self-repellent interac- 
tions. The scattering function is given by 

S{ci)^^j^ dtj^ (exp(*q.(r(t)-r(T)))dT, (16) 

where q is the scattering vector. For a chain with the 
Hamiltonian H{y) and the Green function G(Q), the av- 
erage value ($) of a functional $ reads 



($) = 
(<f>Q = 



$)qG(Q)dQ, 

^ ML)=Q 
G(Q) Jr(0)=O 



$(r(t)) 



(17) 



where I? is a measure on the set of vector-functions y(t ) 
obeying the boundary conditions in the path integral |33| . 
Combination of Eqs. H16|l and (|17|l implies that 



5(q) 



where 
S(q,Q) 



1 

1? 



r(L)=Q 

exp 

r(0)=0 

dt I exp 



S(q,Q)dQ, 

H{v{t)) 



(18) 



(zq • (r(i) - r(T) 



.(19) 

/o Jo ^ ^ 

According to Eq. (|16|l . the scattering function S'(q) sat- 
isfies the normalization condition 



5(0) = 1. 

The radius of gyration of a chain 



(20) 



1 

2L2 



dt / 

"'0 



r(t)-r(r))')dr 



(21) 



is connected with the function S'(q) by the formula 



Rl 



|a5(0), 



(22) 



where A stands for the Laplace operator. For an isotropic 
scattering function that depends on only, S{q) 
S{q^), Eq. is substantially simphficd, 

Rl = -dS'{Q), (23) 



where the prime stands for the derivative. 

IV. EXPANSION OF THE SCATTERING 
FUNCTION 

Expanding the exponent of the Hamiltonian H into the 
Taylor series in we find from Eqs. Q and 115() that 



n— 



H^{r)\^{-ir (vb\ 







X 



'^tj^ p(t-T)(5(r(i)-r(r))dT 



(24) 



Inserting expression (|24|l into Eq. H19|l . changing the 
order of integration and summation, and using Eq. I|17|) . 
we obtain 



^«..Q)^^E^(1F)" 

J dt J exp^iq • ^r(t) — r(T) 



dt / pit-T)S[vit)-r{T)\dT } ,(25) 

Jo ^ ^ / Q 

where (•)q stands for averaging over configurations of a 
chain with Hamiltonian |^ . To derive an explicit expres- 
sion for the nth moment (n = 0, 1, . . .) 



Mn = 



L nL 

dt / cxpf • ( r(i) — r(T 



Jo 

L i-L 



dt / ,5(r(t) -r(r))dr 

^0 



(26) 



Q 



we use the identity 



dt 
Jo 

dti / dri 

^0 



p{t-T)5(v{t)~v{T)) 



dr 



L 

dtn / dr„ 
Jo 



X n ~ T-„i)(5(^r(t„i) - r(r„)^ , 

m— 1 

insert Eq. ((3J, and obtain 

X / dtn / dT„ / dtn+1 / dT„+i TT p{t,n - T„i) 

Jo Jo Jo Jo 



ri+1 

X { exp ^-i ^ k,„ • ( r{tm) - r{Tm) 

m—1 



m— 1 


Q 



= where 



tn+l — t, Tn+l — T , k„+i — — q. 



(27) 
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Ordering the set of points {tm,Tm} as ti > ti > t2 > T2 . . . > in+i > Tn+i and bearing in mind that for any 
integrable function S{ti,Ti, . . . , t„+i), 



dti / dri 

Jo 



Atn+i / S'dr„+i = (2n + 2)! / dti I dn . . . / dt„+i / S'drn+i, 



we find that 



M„ = 



(2n + 2)! 



m— 1 

n+1 



n / dk„, / dti / dr: 



(27r)"^ 

X {^exp|^-i ^ k„i • {v{t,n) - r(T„j) 

m— 1 

For an ordered set of points {tmiTjn}, the average in Eq. I|28|) reads 

ri+l 

fexp|^-i ^ k„ • (^r(t„) - r(T, 



JO 



dT„+l W p{t„i - T„i) 



(28) 



Q 



= n^^^pf-^^' • ^(T) • ki - }Q/k, • t), 



(29) 



where Qi is the Zth component of the vector Q in a Cartesian coordinate frame, kmi is the Zth component of the vector 
km, and T„i — tm — t„i {m — 1, . . . ,n + 1). The {n + l)-vcctors k; and T and the (n + 1) x {n + l)-matrix C(T) are 
given by 



kii 
hi 



k; 

-qt 



c = 



knl 
kn+1 I 

Ti{l-Ti/L) -^TiT2/L 

-T1T2/L T2{1~T2/L) 

-TiT^/L -T2T3/L 



Ti 
T2 

Tn 
Tn+1 



T 

t - T 



(30) 



— TiTn/L —TiTn+l/L 

-T2XnlL -~T2T„+i/L 



-TiT^/L -T2TJL 

— TiTn+l/L —T2Tn+i/L 

C -T{t-T)/L 

-T{t-T)/L {t-T){l-{t-T)/L) 



—TnTn+i/L r„+i(i — r„+i/L) 



(31) 



Taking into account that 

d 

I cxn I — 

2d 



[]cxp(-|^k, . C(t) . k, - jQM ■ t) = expL ^fl^k, . C(T) . + 1(Q; - zg;^(t _ r))k, • T 



1=1 



X exp 



,2d 
1=1 

t — T f boLq^ 



L \ 2d 



t - T 



iQ, q 



(32) 



we find from Eqs. (|2Hl and that 



(2n + 2)! 
1^ 



\nd 



1^1 



n dk, 



q ; All I dri 
/o Ja 



dt I drcxp 
Jo 



771 — 1 



- T„i) exp 
I 



t — T f boLq 



L V 2d 



1 - 



t-T 



L 



iQ q 



where we replaced integration over k„i with integration culating the Gaussian integral over kj and bearing in 
over kj and returned to the initial notation t and t. Cal- mind that /J"" dt /J (p{t - r)dr = /(["(rn - z)tf{z)dz for 
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any function ip{t — t), we obtain 



dti / dri . . . / dtr, I dr„. 



X / (r„ - z)i?„(T,z)dz, 
/o 



(33) 



where 



X exp 



_!^q;,)'T.C-i(T).T 



h r ^ 



and we omit the arguments q and Q of the function R„ 
At n = 0, Eqs. (ESI and ^ read 



y dt J exp^iq • ^r(t) — r(T) 

/ (L-z)i?o(T,^)dz 
Jo 



Q 



with 
i?o(T,z) = 2 exp 



z /boLq^ 



LK 2d V^-l'-^Q-^ 



(35) 



(36) 



To reduce the number of integrations in Eq. H33|) . we 
infer by induction that for any ri = 1, 2, . . ., 



dti / dri 

"'0 



xi?„(T,z)dz = 



dtn / dT„ / (t„ - z) 
"'0 "'0 

L /-L-Ti 



1 



dTi 



dTs 



(n + l)!yo Jo 

X / dr„ / " 



2 ■ 











m— 1 



\ n+1 

-z i?„(T,z)dz. 



Introducing the new variables 



(37) 



L-^Tfc (m= l,2,...,n), (38) 



we present Eq. I|37|l in the form 



dti / dri 
Jo 



xi?„(T,z)dz : 



dtn / dT„ / (t„ - z) 

Jo Jo 

, dui / du2... 



1 



nUn—1 PUn 

X du„ / (M„-z)"+ii?„(u,z)dz, (39) 

Jo Jo 

where i?n(u, z) is obtained from i?„(T, z) after substitu- 
tion of the expressions 



Ti = L-ui, Tm^u,n-i~u,n (m = 2,...,n), (40) 

corresponding to the inverse transformation. To find an 
exphcit expression for the function i?„(u, z), we use the 
formulas 



detC(T) = (l - ;^ Tm) n 

m—1 m—l 

1 " 

= — (L - ui)'u„ (u„_i - u„), 

m=2 

n n 

T • C-i(T) • T = L ^ - ^ T„ 

m — l m—l 



(41) 



that follow from Eqs. (j^OI), ((HU and gOll. Insertion of 
expressions (|41f) into Eq. H34|l results in 



i?„(T,z) = (2n + 2)!(^) ^"i^ 

p(L - Ui)p{ui ~U2) ... p{Un-l - Un) 
X ^ 

(L - Ui)2{ui - 7/2)2 . . . (li„_i - Un)"^Un 

d{L - u„) 



X exp 



2boLun 



■(q 



z6o ^2 



z /boLq^ 



(42) 



It follows from Eqs. (|23l, CHI), (El, (|S3l, JSEI), (EHJ and 
(E^ that 



S(q,Q) 



Go(Q) 



2 / {L — z) exp 



L 



1 - 



z 



27r6o 



p{L — ui)dui 
(L — ui)^ 



2d 

p(wi - U2)du: 



- ) - • q 



dz + 



(-1)" 



n=l 

p(u„-i - u„)du„ 



X exp 



(i(L - U n) ^ 
2boLu 



(ui-U2)2 Jo (u„_i -u„)5u^! -^0 

2 r /&oLg2 



•(2n + 2)! 
2i y (n + 1)! 

(w„ - 2)"+' 



*n)/„ ibo \2 z (hoLq^ ( z\ 



dz 



(43) 
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V. THE LAPLACE TRANSFORM 

The Green function for a Gaussian chain in a d- 
dimensional space reads 



Equation (|44() differs from Eq. ^ by the pre-factor G. 
This coefficient will be found later from normalization 
condition H20|l . Substitution of expressions H43|l and H44|l 
into Eq. (|18|) results in 



2dL 



dz / exp 



2boL 



d-l , 



(2n + 2)! 

{n + 



2)!/ d 

1)! V2^6o/ Jo 



p{L — ui)dui 



p{ui - U2)du2 



X cxp 



6og^z(u„ - z) 

2dUn 



{L - ui) 2 Jo (mi - U2) 2 
dz / exp -— \ Q-q dQ 

V 20oU„ Un ' 



Calculation of the Gaussian integrals yields 



5(q) 
G 



D 



boL 2 



ra=l 



J2! 

2i 



(-1)" fvb^-\''^i2n + 2)\ 
in + l) 



1)! V27rfeo/ 



— ui)dui 



(i - Jo [ui - U2) 



p{ui - U2)du2 



f 

Jo 



(u„ - 0)"+-^ exp 



r 



where 
Sd 



boL 2 



2(i 



1 = cxp 



exp 



2d 

boW_ 
2d 



zj {L — z)dz 
(l-a;))a;da; (46) 



is the Debye function. Introducing the notation 
P(i,q)=i2 



Sin) ^ f^oL 2\ 



we infer that 



G 



(-1)" /i'6^"\"(2n + 2)! 



(47) 



p{L — ui)dui 



(n + 1)! \2TTbo 



p{ui - U2)du2 



{L-uiY 



[ui - ^2)2 



(u„_i - U„)^ 



(u„ - 0) 



n+1 



X exp 



/_6oV\ 
V 2d ) 



dz. 



(48) 



2d 



dz. 



(45) 



Denote by P(A, q) the Laplace transform of the func- 
tion P(I/,q) with respect to L. Multiplying Eq. by 
exp(— AL) and integrating over L from zero to infinity, 
we find that 



n—1 

00 

exp(— AL)dL 

p{ui - U2)du2 
{ui - U2)^ 



(-1)" p6g-S»(2n + 2)! 
2L / (n + l) 



1)! V27r6o/ 



p(L — ui)dui 
[L — ui)i 



u„)dw„ 



(u„_l - Wn)^ 



dz. 



(49) 



According to the convolution theorem for the Laplace 
transform, for any function (^(u). 



exp(— AL)d_L 



p{L - ui) 



{L - ui) 
A(A) / <y9(iti) exp(— Aiti)dui 



-Lp{ui)dui 



To derive an explicit formula for the function P{L,q), 

the following technique is employed: (i) first, we perform where 

the Laplace transform of Eq. I|48(l . (ii) afterwards, we 

sum the infinite series in the transformed space, and (iii) 

finally, we perform the inverse Laplace transform. 



A(A) = / p{x)x 2 exp(-Aa;)da 



(50) 



(51) 
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Applying Eq. (|50|l n times, we obtain 

exp(-AL)dL PiL^^ 

p{ui - U2)du2 /-""-i p(un-l - Un)dUn 



X 



(Ul - U2)2 

A"(A)/. bnQ^\-^ 



(U„_i - Un)^- 



Combination of Eqs. H49|l and (|52|l implies that 



P(A,q) = A-^ A + 



2d 



f (53) 



with 



V /dbo 

The pre-factor in Eq. H53() reads 



(54) 



(2» + 2)! , ,, r(2(n + l)) 

fln = i = 2 n + 1 — — — — . 

n\ 1 (n + 1) 

Application of the Legendre duplication formula for the 
gamma-function yields 

4"+! 













-1) 




1)1 











Insertion of this expression into Eq. (|53ll results in 

- oo 



1 oo „ • 

x{2vr ~ -j2{~irT(n+ -){2vr 

n=l 

where we use the shortcut notation 



A 



(55) 



Bearing in mind that for any X 

^ oo „ 



2 L 4 



we find that 



3^/^ 



0FA2 

+(i + ^)E(-irr(n + ^)(2y)" 



(56) 



With reference to the generalized method of Borel's sum- 
mation W5l , an infinite series 



is associated with the function 



a{z) = cxp(-z)E 



•nA- — 



rir(n + |) 



(57) 



(58) 



Series H57f) is Borel summable, and its sum equal uq, pro- 
vided that the integral /p a{z)dz exists for any s > 
and converges to ctq as s ^ 00 in the classical sense. 
The summation method with function (|58|l is regular: 
for any convergent series H57|) . ctq = Yl^=i ^n- Applying 
the above technique to the infinite series in Eq. (|56|l . we 
find that for any X , 



E(-irr(n+-)x^ 



exp(-z) /^"^ dz, 
1 -f- ^ z 



where we performed summation of the geometrical series. 
The expression in square brackets in Eq. (|56|l reads 



V 



= V 



-2(1 + V) 
1 - 2z 



exp(— z)z2 



/o 1 + 2Vz 
where we used the identity 
'■^ 



1 + 2V"z 
exp(— z)z2 dz, 



dz 



(59) 



exp(— z)z2 dz. 



Inserting expression 159|l into Eq. H56|) and using Eq. 
(|55|l . we arrive at the formula 



^^(A,q) 



_4uo_ / 

./^aH 



A + 



2n -1 



0FA2V' ' 2d 

A(A)(1 - 2z) 
A + 2woA(A)z 



exp(— z)z2dz (60) 



for the Laplace transform of the function P(L,q). Our 
aim now is to find the function 



P(i,q)=/:-MP(A,q) 



where £ ^ stands for the inverse Laplace transform. Tak- 
ing into account that 



(A. 



2d 



2n -1 



= cxp 



2d 



and applying the convolution theorem for the inverse 
Laplace transform, we obtain 

A /"OO 

P(L, q) = / exp(-z)(l - 2z)zidz 
\^ Jo 

boiL - l)q^ 



X ; exp 
'0 



\ 2d ) 



C-^{w{\,z)^^ dl, (61) 
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where 



W{\z) 



A(A)A- 



A + 2t;ozA(A)' 



(62) 



As the inverse Laplace transform in Eq. Ht)l|) cannot be 
determined in the general form for an arbitrary integer 
d, we consider particular cases d — 1 and d = 3. 



VI. A SELF- AVOIDING CHAIN IN AN 
ONE-DIMENSIONAL SPACE 

As the integral in Eq. H51|) with d=l converges with- 
out regularization, we set p — \ and find that 



A(A) = j X 2 exp(— A.T)da:; = -^J^. 
Insertion of expression into Eq. results in 



(63) 



A2 + (j){z) 



(z) = 2vq^/ttz ~ vz 



(64) 



Taking into account that for any positive a, (3, and 0, 



where 



n=0 + P) 



is the Mittag-Lefller function we infer that 
Combination of Eqs. H61|l and (|67|l yields 



P{L, q) = 4wo / exp(-z)(l - 2z)zidz 
Jo 



(65) 



(66) 



(67) 



(_M^)^3,(-,(.).i).i 



dL 



Setting X = l/L and substituting expressions H54|) and 
()64|1 into this equality, we arrive at the formula 

P(i,q) = -^/ exp(-z)(l -2z)z5dz 
V"" Jo 

^-^^(l^x))Es^.i-Xxh)xUx, (68) 



X / exp 
/o 



where 



2L 
V bo 



(69) 



Insertion of Eq. (|^ into Eq. (gYJ yields 

-[ 2y I b'^q'^ 

^(q) = G[5.(^)+-|y^exp(-^(l-.) 

/>oo 

xa;^dx / exp(— z)_E3 i (— z)(l — 2z)z2 dz 
Jo ^ ' ^ 

It follows from Eq. that for any c > 0, 



(70) 



exp(— z)£^3 7(— cz)(l — 2z)z2dz 



t^o r(i- + i) 



(71) 



Substituting expression with d = 1 into Eq. l(7n|) 
and using Eq. we obtain 



'^(q) = exp(^-^^(l - a;))xdx 



■n— ^ ^ 



(72) 



Setting q = in Eq. H72|) and employing normalization 
condition (|20|l . we infer that 



G = 



- r(n+§)(n + l) 

r(|n + 3) ^ • 



(73) 



Taking into account that for any a; > 0, 

- r(n+|)(n+l) 
„t^o r(|n + 2) (-^^ 

/tT /""^ 3 

"2/0 ^t2(-^y)c^p(-y)ydy, (74) 

we derive an explicit presentation of the scattering func- 
tion in terms of integrals of the Mittag-Leffler function 



exp( 7^i^ ~ x) \xdx 



ir 2(~^y) exp(-t/)ydy 



(75) 



To determine the radius of gyration i?g, we insert Eq. 
(|72|l into Eq. (|22|l . calculate the integrals and find that 



2 262(5^ r(n +§)(«+ 1) 



Ri = 



By analogy with Eq. (|74|l , one can write 



(-X)"- (76) 



t^o r(i" + 4) 

2 



= ^ I Ei^{-xy)eM-y)y'^y- (77) 
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0.0 



8.0 



FIG. 1: The dimensionless ratio ^ versus the dimensionless 
strength of segment interactions x- Solid line: exact solution. 
Circles: linear approximation. 



Combination of expressions (|73|l , 174|l , (|76|l and (|77(l im- 
plies that 



El A-XV) exp(-2/)2/d?/ 



El J-xy) exp(-y)ydy 



(78) 



At weak self-avoiding interactions, when the dimension- 
less strength of segment interactions x is small compared 
with unity, we disregard terms beyond the first order of 
smallness in Eqs. (173(1 and ((77|l to obtain 



32x 
6 V" ' 105V^ 



(79) 



Equation H79I) shows that at any x > 0, the radius 
of gyration of a weakly self-repellent chain exceeds that 
of a standard Gaussian chain, Rgo — \h^- To evaluate 
the effect of excluded-volume interactions on the radius 
of gyration, we perform summation of the infinite series 
(|73|l and H7t)|l numerically. The dimensionless ratio ^ = 
Rl/Rl is plotted versus the dimensionless strength of 
segment interactions x in Figure 1 , which shows that the 
radius of gyration increases monotonically with strength 
of self-repellent interactions. Linear approximation H79|l 
ensures a correct fit of the dependence Rg{x) at relatively 
small strengths (x < l)i but overly estimates the radius 
of gyration at higher values of x- 



VII. A SELF-AVOIDING CHAIN IN A 
THREE-DIMENSIONAL SPACE 

For a self-repellent chain in a three-dimensional space 
(d = 3), the use of a regularizing function p{t) is neces- 
sary, because the integral with p = I in Eq. H51|) diverges. 



It follows from Eqs. (O and lfET|) that 

A(A) = 2^/¥(VAT7 - \/a) , 

where r = L^^. This relation together with Eq. Ht)2|) 
results in 



WAX,z) 



VA + r - \/A 



X + ^{z){VxTV - Vx)' 

Bearing in mind that for any X, 
X 



(80) 



^(-0)"-iX", 



we infer that 

oo 

iyi(A,z) = 5](-</.)"-i( 



l + <j>X ^ 



A 



(81) 



To find the inverse Laplace transform of the function 
W{X, z), we begin with the formula 



V / s — >x 



(82) 



where 



^»W-E ^,r(L + n + l) (2) 



is the Bessel function of the first kind of order n, and 
a > is an arbitrary constant. We now apply the well- 
known assertion which says that if <I>(s) is the Laplace 
transform of a function (l){x), then 



dz. 



Applying this proposition to Eq. I(82|l . we obtain 



na 2 



J, 



.(^/az)exp(-—)z Mz. (84) 
V 4a;/ 



Inserting expression ((83|l into Eq. H84|l . performing in- 
tegration term by term, and replacing the resulting se- 
ries with an appropriate Mittag-Leffler function, we infer 
that 



n 2^ 

na X 2 ^/n 



20F 



2 / Jl\ 21 

— r(-)i?i^_„+i(-ax). (85) 
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Bearing in mind that 



we arrive at the formula 



£-1 



S'^ 2 / S — >X 



r(n-i)' 



and applying the convolution theorem for the inverse 
Laplace transform, we find that 



C 



'\Wi{\,z) 



X / a;"(l — x)™ ^dx 



X / y Ml - 2/) ^ exp{-rlxy)dy. 
Jo 



(87) 



1 



+1 r(n — i) 



X r(.-0"-l/:-^((^-^)'^ 



dt. 



Substitution of expression (|85() into this relation yields 



r(t) 



X / {x-tr-H^El.^^,{-at)At. (86) 
Jo 



Combination of Eqs. H81|l and (|86|l results in 
C-^(w,{\z)) 



X / x— (1 -x)""^£;f„+i(-Wa;)dx, 
/o 



where we set x — t/l. Using the integral presentation of 
the Mittag-LefFler function 



1 



r(f)r(| + i) 



X / ^(1 - y)S exp(-r/a;y)dy, 

Jo 



Application of the convolution theorem for the inverse 
Laplace transform to Eq. H8UI) results in 



-'{W{\,z)) 



= 2^ j {l-t)C-^{Wi{\z)^ ^ At, 

where we employed the formula C~^{X~'^)\^x = x. Set- 
ting u = t/l\n Eq. (|55|l and substituting expression l|S7|l. 
we find that 



(to+ 1)(-W50)™ 

^^r(m + i)r(f + 1) 



X / u^"T~(l — w)dM / x^ (l ~ x)™^ dx 

Jo Jo 
pi ^^^^ 
X / y~2— (1 _ exp(— r/-ua;y)dy. (89) 

Jo 



Combination of Eqs. l|?)T|) and yields 



ivor {m + l){—r) 



E 



/^i^or("^ + i)r(T + i)^o 



exp(-z)(l-2z)zi0"(z)dzy exp(- ) 



l^-dl 



X j u ^ (1 — u)du / a; 2 (1 — a;)™ ^dx y 2 (1 _ y) 2^ exp(— rZMxy)dy. 
Jo Jo Jo 

Introducing the notation t = l/L and inserting expression (|54|l and (|80|) into this relation, we find that 



2xL 



E 



i-xY 



, exp(-z)(l-2z)z™+5dz / expf- — (l-t)V^dt 



X / u 2 (1 — u)du / X 2 (1 — x)™ 2da; / y 2 (1 — y) 2 exp{—K,tuxy)dy, (90) 
Jo Jo Jo 

I 



where 



(91) 



The integral over z is calculated explicitly by means of 
the gamma-function. 



exp(-z)(l - 2z)z"+3dz 



'2{m + 2)T(^m+ - 



(92) 
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Inserting expressions (|46|l and H90|) into Eq. (|47() and 
using Eq. I|92|) and elementary properties of the gamma- 
function, we arrive at the formula for the scattering func- 
tion of a self-avoiding chain in a three-dimensional space 

S{q)^2G f exp(-^(l 



E 

m=0 
1 



(to + i)(TO + |)(to -I- 2)(-xi5)'"+i 



r(f + i) 

1 



3m+l 



X I u ^ {1 — u)du / a; 2 (1 — x)™ 
»i 

1+1 



X / y 

/o 



(1 — y) 2 exp(— KiMxy)dy 



(93) 



The function (|93|l is determined by three parameters: (i) 
the mean-square-root end-to-end distance of an appropri- 
ate Gaussian chain 6, (ii) the dimensionless strength of 
segment interactions a-nd (iii) the dimensionless length 
of a chain k. The pre-factor G is found from normaliza- 
tion condition H20() . 



A. Weak excluded-volume interactions 

At weak self-avoiding interactions between segments 
(x ^ 1), we disregard all terms beyond the first order of 
smallness in Eq. H93|) and obtain 



2G I exp(^ |-(l-t)jMt 

1 ^ pi ^ 

1 — "'^'^^ I 1/2(1— u)du / (1 — x)^2clx 



X 



71" V' Jo 







X / y 2 (1 — y)2 exp(—Ktuxy)dy 



(94) 



It follows from Eq. that 



y 2 (1 _ y)2 exp{-Ktuxy)dy = -El^i-i^tux). 
Taking into account that for an arbitrary positive z, 







we find from Eq. that 

5(q) = 2G'/ exp(-^(l-<))M< 



E 



1 ' 6V 
^ 6 
r(n+ i)(-«;t)^ 



i (1 — ii)dM 



2V^^or(^ + 2)r(n+|)7o 
Calculating the internal integral, we arrive at the formula 



5(q) = 2G / exp 



^ 6 



62 g2 



(1-t) 



E 



r(n+ i)(-Kt)" 



2V^,'^or(" + 2)r(r 



Mi. (95) 



To determine the coefficient G, we set q = in Eq. H95|l , 
substitute the obtained expression into Eq. I|2U|) , perform 
integration, and obtain 

r(n + i)(-K)" -"-^ 



G = 



o — ^ 

1-^E 



7r;^^r(n + 2)r(n+f) 



It follows from Eqs. and ((23 that 



6 



1 



9xg r(n-fi)(-«)^ 



^/^^g^(n + 2)^(n+^) 



(96) 



(97) 



where the integrals are calculated term by term. Com- 
bining Eqs. (|95|l and (|96|l and neglecting terms beyond 
the first order of smallness with respect to x, we find that 



6 



3x y^ r(n+l)(n+f)(-^)- 



n=0 



(98) 



Bearing in mind that 

r(n+i)(n + f) 



r(n-h2)r(n+ ^) 



and 



^I/3(-«2;)da; 



1 



E 



r(n + f) + r(n + 1) 

(n + l)r(n+^)n! 

r(n + 7) (-^)' 



r(7) (" + i)r(«« + /3) 



we present Eq. 197|l in the form 



3x 



(99) 

Equation (|99|) provides an explicit expression for the ra- 
dius of gyration of a self-repellent chain in the limit of 
small strengths of segment interactions (x 1) for an 
arbitrary dimensionless parameter k. 



B. Asymptotic expression for the scattering 
function 



Our aim now is to derive a formula for the scattering 
function at strong self-repellent interactions with w — > oo 
under the assumption that k — > and the product vn in 
Eq. (|91|l remains finite. 



3uK 
lim 

V — >oo, K — >0 TT 




2x>0. 



(100) 



The coefficient "2" is introduced in Eq. H100|) for conve- 
nience of further transformations. 

Replacing the exponential function in Eq. H93|l by 
unity and calculating the integrals in terms of the 
gamma-function, we obtain 

S{ci)^^ j\x^(~^^{l-t))tdt 

^{n + l)r(in + \)T(\n + l)r(n + |)(-2xii)" 



n=0 



r(n+ l)r(|n + 2) 
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where we included the first term in the square brackets 
into the infinite sum. Applying the Legendre duplication 
formula, we find from this relation that 

r|-r(„ + i)(,. + i) 1 

r |n + 2 V A y J V / 

'-ri=0 ^ 1 ' -I 

Comparison of Eqs. H72|l and (|101|l implies that in the 
limit (|100|) . the scattering function for a self- repellent 
chain in a three-dimensional space coincides (after the 
natural replacement iq^ — > \q^^ with that for a self- 
avoiding chain in an one-dimensional space. Based on 
this similarity and using Eq. (|75|l . we arrive at the closed- 
form expression for the scattering function 

„\ , 
5(q) = exp(--52g2(i_t))Mt 

/o°° ^\ S-X^^y) exp{-y)ydy 

(102) 

that extends the Dcbye formula for a standard Gaussian 
chain to self-avoiding chains. It follows from Eqs. H23|) 
and (|102|l that the radius of gyration of a self-repellent 
chain in a three-dimensional space is given by Eq. (|78|l . 
and its dependence on strength of self-avoiding interac- 
tions is presented in Figure 1. 
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FIG. 2: The Guinier plots of the scattering function S for 
a standard Gaussian chain (x ~ 0.0) and for self-repellent 
chains with various strengths of excluded-volume interactions. 
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FIG. 3: The scattering intensity / versus the amplitude of 
scatte ring vector q. Symbols: experimental data on a PDMS 
blend Q at volume fractions of H-PDMS 0.06 (unfilled cir- 
cles), 0.26 (filled circles), and 0.52 (asterisks). Solid lines: 
results of numerical simulation. 




FIG. 4: The scattering intensity / versus the amplitude of 
scattering vector q. Large symbols: experimental data on 
solutions of PS/PBLG [isll in dioxane (unfilled circles) and 
dioxane/TFA mixture (filled circles). Solid lines: results of 
numerical simulation. Small circles: approximation by the 
Debay function with b = 137.2 nm. 



VIII. RESULTS OF NUMERICAL SIMULATION 

To perform numerical simulation, we return to presen- 
tation (|101|l of the scattering function and write 

^(q) = (^c„(-x)")" / e^J--b'q^l^t)) 

CO 

xj^ani-xt^rtdt, (103) 

n=0 
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FIG. 5: The scattering intensity I versus the amplitude 
of scattering vector q. Symbols: experimental data on a 
PDME/PEMS blend at the temperatures T = 100, 200 and 
250 °C, from top to bottom, respectively l50|l . Solid lines: 
results of numerical simulation with x — 0.84 . 



FIG. 7: The mean-square-root end-to-end distance for an 
equivalent Gaussian chain b versus temperature T. Symbols: 
treatment of observations on PDMS/PEMS blends. Solid 
lines: approximation of the experimental data by Eq. (104). 
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FIG. 6: The scattering intensity / versus the amplitude 
of scattering vector q. Symbols: experimental data on a 
PDME/PEMS blend at the temperatures T = 150, 200 and 
300 °C, from top to bottom, respectively [sof. Solid lines: 
results of numerical simulation with x = 0.67 . 



where 



r(|n-f 2) 



2a„ 



3n 



We sum numerically series (|103|l for several values of x 
and present the dependence S(q) in Figure 2 (the Guinier 
plot). Given an amplitude of the scattering vector q, an 
increase in strength of self-repellent interactions results in 
a pronounced decrease in the scattering intensity. With 
the growth of x, the graph of the scattering function flat- 
tens at relatively small values of q and it becomes more 
close to the graph of the Giunier function. The latter 



function is believed to adequately describe the scattering 
intensity of chains with strong excluded- volume interac- 
tions [43. 

We proceed with matching observations [small angle 
neutron scattering (SANS)] on poly(dimethylsiloxane) 
(PDMS) blend 5HC/11DC with various volume fractions 
of hydrogenated PDMS (H-PDMS) reported in |4||. A 
description of the experimental procedure is given in '48| . 
Assuming the scattering intensity I(q) to be proportional 
to the scattering function S{q), I = KS, where K is an 
optical constant |i32i] . we approximate the experimental 
data depicted in Figure 3 by using three adjustable pa- 
rameters: b, X s-nd The best-fit quantities b and x 
are determined by the steepest-descent method, and K 
is found by the least-squares algorithm. Figure 3 demon- 
strates an excellent quality of matching the observations 
with X varying in the interval between 2.3 and 4.3 . 

To compare the accuracy of fitting experimental data 
by Eq. H1U3|) and by the Debye function H46|l , we approx- 
imate observations (SANS) on poly (styrene-dg )-b-poly (- 
benzyl L-glutamate) (PS/PBLG) reported in The 
scattering intensity I is plotted versus the amplitude of 
scattering vector q in Figure 4 for solutions of PS/PBLG 
in dioxanc and mixture of dioxane and trifiuoroacetic acid 
(TFA) with 20 wt.-% of TFA. For a detailed descrip- 
tion of the experimental procedure, see where it is 
demonstrated that the presence of TFA in the solution 
results in transition from the rod-coil architecture of the 
diblock polymer to coil-coil topology of macromolecules. 
Figure 4 shows that the best-fit parameter x vanishes for 
the coil-coil system, and Eq. (|1U3|I is reduced to the De- 
bye function. For the rod-coil architecture of the diblock 
polymer, x equals 3.1 (approximately), which means that 
the effect of excluded-volume interactions between rods 
is substantial. The Debye fit of the observations on so- 
lution of PS/PBLG in dioxane is rather poor compared 
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with the approximation by Eq. (|103|l . and the best-fit 
value of b found when self-avoiding interactions are dis- 
regarded appears to be too large. 

To study the effect of temperature T and mass-average 
molecular weight Af^ on the material constants in Eq. 
()103|l . we approximate observations (SANS) on solu- 
tions of symmetric blends of poly(dimethylsiloxane) and 
poly(ethylmethylsiloxane) (PEMS) in toluene. A de- 
tailed description of polymers and the experimental pro- 
cedure is given in |50|. The experimental data and the 
results of numerical simulation are presented in Figures 
5 and 6 for two blends at temperatures ranged from 100 
to 300 °C. To reduce the number of adjustable parame- 
ters, wc find X for each blend by matching observations 
at the lowest temperature and use this quantity with- 
out changes to approximate the scattering intensity mea- 
sured at other temperatures (this means that b is the only 
adjustable parameter in the fitting procedure). Figures 
5 and 6 demonstrate good agreement between the ob- 
servations and the results of numerical analysis (in con- 
trast with the model used in the original work |5^] that 
contains the same number of adjustable constants, but 
shows rather poor quality of matching the experimental 
data). The parameter x equals 0.84 for the blend with 
low molecular weight and 0.67 for that with high Mw. 

The mean-square-root end-to-end distance for an 
equivalent Gaussian chain b determined by fitting these 
data is depicted in Figure 7 together with its approxima- 
tion by the Arrhenius-type dependence 

\nb^ f3o + (3iT-\ (104) 

where the coefficients f3Q and (3i are found by the least- 
square technique. Formula H104|l is employed as a phe- 
nomenological relation that describes the effect of tem- 
perature on the segment length 6o = b^/L. Figure 7 
reveals that (i) the number of statistical segments in- 
creases with molecular weight, and (ii) given M^, the 



Kuhn length decreases with temperature (the growth of 
T is tantamount to an increase in the intensity of thermal 
fluctuations, which, in turn, implies that the characteris- 
tic length of correlations between monomers diminishes). 



IX. CONCLUSIONS 

Introducing regularizing function p in the form jSJl, 
we have derived explicit expressions for the scattering 
function S'(q) of a self-avoiding chain in a d-dimensional 
space. Analytical relations are developed for the form 
factor and the radius of gyration of a self-repellent chain 
in the one-dimensional and three-dimensional cases. Re- 
sults of numerical simulation demonstrate that the con- 
ventional linear approximation for the radius of gyration 
(weak excluded-volume interactions between segments) 
ensures an acceptable agreement with the exact solu- 
tion when the dimensionless strength of interactions is 
less than unity. Fitting of experimental data (small- 
angle neutron scattering) on several polymer solutions 
reveals that the values of the dimensionless strength of 
self-repellent interactions exceeds the ultimate value at 
which this approximation is reliable. 

A novel approximation for the scattering function is 
proposed (high strength of segment interactions v oo 
and large characteristic length of internal inhomogene- 
ity Lq/L ^ oo). In the limit H100|) . the scattering func- 
tion for a self-repellent chain in a three-dimensional space 
is similar to that for a self-avoiding chain in an one- 
dimensional space. Equation (|102|1 deduced in this ap- 
proximation provides a simple generalization of the De- 
bye formula H46|l for the form factor of a standard Gaus- 
sian chain. Comparison of numerical results with obser- 
vations reveals that Eq. H103|l adequately matches exper- 
imental data, and its material parameters are affected by 
external conditions in a physically plausible way. 
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